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Sideband cooling is a technique that potentially allows mechanical resonators to be prepared in 
their ground states, important for future applications in quantum technologies. Tian has recently 
shown that side-band cooling can be implemented by modulating the coupling between a nano- 
resonator and a superconducting oscillator, a process of frequency conversion [L. Tian, PRB 79, 
193407 (2009)]. While side-band cooling is usually treated in the steady-state regime, the effective 
resonant coupling will also generate near perfect state-swapping from the superconductor to the 
mechanical resonator. We perform numerical simulations of this system, examining the ground- 
state cooling achieved by the state-swapping. Further, we show that the superconducting oscillator 
can be used to control the amplitude and phase of the resonator, while simultaneously cooling it, 
and thus act as a coherent "quantum feedback controller" . 
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Nano-mechanical resonators can now be built with 
high frequencies and high quality factors [THS] and in- 
terfaced with mesoscopic superconducting circuits [H [5J . 
In order to observe and exploit the coherent quantum 
properties of these devices, one must first be able to 
cool them to their ground state. Ground-state cooling 
is achieved when the average number of phonons in the 
resonator, (n) , (also called the "average occupation num- 
ber" ) is much less than unity. Feedback control via mea- 
surements could in theory be used to perform such cool- 
ing [BJ, but this approach is not yet feasible with current 
technology. Two classes of methods have been devised 
to date to achieve ground-state cooling without using 
measurements. The first involves periodic coupling to 
a superconducting qubit or qutrit [TJ [5] ■ The second con- 
sists of various ways to realize side-band cooling [51lll|. 
also known as "radiation pressure" cooling or "dynamical 
back-action" cooling. Recently Tian [TT] proposed per- 
forming sideband cooling by using a linear coupling be- 
tween a nano-resonator and a "auxiliary" superconduct- 
ing L-C oscillator, and modulating this coupling at the 
difference frequency between the two resonators. This 
formulation is especially nice, because it shows that the 
cooling comes from the fact that the frequency conver- 
sion makes the auxiliary oscillator look to the mechanical 
resonator as if it is much colder, even though it is really 
at the same ambient (background) temperature. 

In a nutshell, sideband cooling viewed as frequency 
conversion can be described as follows. If one couples a 
low-frequency resonator (which we will call the target) 
to a high-frequency resonator (the auxiliary) via a lin- 
ear coupling modulated at the resonators' difference fre- 
quency, then the two resonators appear to each other as 
if they are on-resonance. Since both resonators are at the 
ambient temperature, the auxiliary has a lower occupa- 
tion number (is closer to its ground state) , and therefore 
looks to the target as if it is at a lower temperature. 
If we now ensure that the auxiliary has a much faster 



damping rate than the target, energy that flows from the 
target to the auxiliary is quickly lost to the bath. The 
result is that the target is cooled close to the effective 
temperature of the auxiliary. This is the cooling method 
in a nutshell, minus a couple of important additional de- 
tails that we will discuss below. The above scenario is 
perfectly suited for cooling a nano-resonator by coupling 
it to a superconducting L-C oscillator [TT] or "stripline" 
resonator [TJ [4j [5] . This is because these superconducting 
resonators have both much higher frequencies and much 
higher damping rates than their mechanical cousins. 

The caveat to the above description is that the mod- 
ulated coupling only generates perfect transfer of energy 
quanta between the two resonators to the extent that the 
rotating wave-approximation (RWA) is valid. This is the 
case so long as the coupling rate between the oscillators, 
as well as their damping rates, are small compared to 
the frequencies of both oscillators. If the rotating wave 
approximation is not valid, then the coupling also gener- 
ates quanta in both oscillators, leading to heating. This 
is why side-band cooling only provides good ground-state 
preparation in the resolved-sideband limit, in which the 
damping rate of the auxiliary is much smaller than the 
frequency of the mechanical resonator. 

In addition to providing sideband cooling, a resonant 
linear coupling between two oscillators has the remark- 
able property that it generates a perfect state-swap be- 
tween the oscillators for an evolution time of t = ir/ (2g), 
where g is the interaction rate (defined precisely be- 
low) |23j . Since the auxiliary oscillator is initially in its 
ground state (at the ambient temperature) , a single state- 
swap will prepare the target in its ground state, thus 
providing a second way to achieve ground-state cooling. 
We examine both methods of cooling here, via numerical 
simulations. Note that this technique of cooling by swap- 
ping the oscillator states is closely related to the cooling 
method using an auxiliary qubit suggested by Zhang et 
al. [TJ. The present method is much more efficient, how- 
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ever, as the auxiliary oscillator can absorb all the energy 
in a single interaction time t, while on this time-scale an 
auxiliary qubit only absorbs a single quanta. 

In the second part of our analysis, we show that the 
frequency conversion coupling can be used to control the 
coherent state of the resonator, while simultaneously pro- 
viding sideband cooling. In this way, the L-C oscillator 
reduces the noise in the target system as well as con- 
trolling it in real time. This is precisely the action of a 
feedback control loop — note that the noise that is elim- 
inated by the cooling process does not have to be ther- 
mal in origin, it is equally effective at extracting noise 
from any source. In the usual feedback loop, measure- 
ments are made on the system, and actions taken based 
on these measurements to reduce the noise and effect 
control. The frequency-conversion coupling with the L-C 
oscillator is therefore an example of a coherent feedback 
control loop [T6H20] . in which the information flows di- 
rectly from the target system, to another quantum sys- 
tem that acts as a controller, without the use of mea- 
surements. There is therefore a fundamental connection 
between sideband cooling, state-swapping, and coherent 
quantum feedback control. 

Expressions for the average occupation number achiev- 
able by sideband cooling in the steady-state have been 
derived in [TTJ [T2] . We note that one can also derive a 
very simple approximate expression for the steady-state 
cooling, and we present this now for comparison with 
the numerical results presented below. In what follows 
we will denote the frequencies of the target and auxil- 
iary systems as u> and f2, respectively, and their damping 
rates as 7 and n. We will also denote the coupling rate 
between the two systems as g, and assume that the av- 
erage occupation number of the auxiliary is negligible at 
the ambient temperature (this is true for superconduct- 
ing resonators with fi/(27r) > 5 GHz in dilution refriger- 
ators). If we take the coupling constant g to be at least 
as large as the auxiliary damping rate k, then it seems 
reasonable that the energy flow rate out of the target will 
be limited by, and thus approximately equal to, k. In this 
case equating the heating and cooling rates, one obtains 
an estimate for the steady-state occupation number: 

(n) c « [7/(7 + «)]"t, (1) 

where tit is the occupation number at the ambient tem- 
perature T. An example for realistic parameters would 
be a 20 MHz nano-resonator with quality factor of 10 5 
(7/277 = 200 Hz) coupled to an L-C oscillator with a fre- 
quency of 5 GHz and a damping rate of K = 10 6 s _1 [2"2l . 
Choosing a coupling rate of g = 2 x 10 6 s _1 gives 
(n) c « 1 x 10- 3 n T . 

Numerical results for cooling: We now treat the 
frequency conversion, steady-state cooling, and state- 
swapping in more detail. We perform numerical sim- 
ulations to check Eq.Q, and the limiting effect of the 
coupling. We consider the following linear coupling be- 
tween a nano-resonator and a superconducting resonator: 

H = hwo) a + h\(a + a*)(b + tf) + Mltfb, (2) 



in which a and b are the annihilation operators for the 
nano-resonator and stripline, respectively. Here A is the 
strength of the coupling between the two oscillators. Tian 
has shown explicitly how to obtain a linear coupling be- 
tween a nano-resonator and an L-C resonator (TJJ. Al- 
ternatively a mechanical resonator can be coupled to a 
superconducting "stripline" resonator via the nonlinear 
interaction H lnt oc b^b(a + a'), or to an optical cavity 
via the same interaction (in this case it is referred to as 
"radiation-pressure"). By driving the stripline one can 
obtain from this an effective linear interaction, and this 
is the usual context in which radiation-pressure cooling 
has been analyzed in the past. 

To achieve the frequency conversion we now modulate 
the coupling strength at the difference frequency A = 
fl — cd, so that A = gcos(At). Writing the Hamiltonian 
in the interaction picture, the result is 

H in t = ng{atf + a f 6) + hg{abe- l2ut + aW 2 "*), (3) 

where we have dropped all terms oscillating at f2, since 
this frequency is much higher than all the other dynam- 
ical timescales. 

The time-independent part of the interaction Hamil- 
tonian in Eq.([3]) describes the resonant exchange of en- 
ergy quanta between the two oscillators. This is what we 
need for cooling, state-swapping, and coherent control. 
The time-dependent terms will interfere with all these 
processes. In the usual picture of side-band cooling, the 
time-dependent terms are viewed as generating a heat- 
ing rate [TTJ [T^] . But the frequency conversion picture 
we have here does not afford us with a simple interpre- 
tation. The rotating wave approximation (RWA), which 
allows us to drop the time-dependent terms, requires that 
g <C oj. In this case the coupling between the two oscil- 
lators is perturbative, and, in particular, has a minimal 
effect on the ground state of the target. Since g deter- 
mines the rate at which energy is transferred out of the 
target, this ultimately places a limit on the rate of the 
cooling. Nevertheless, since nano-resonators have frenqe- 
cies in the 10's to 100's of MHz, cooling rates can still be 
on the order of 10 6 s _1 . 

We model the thermal damping of both oscillators us- 
ing the quantum-optical master equation. This is an ac- 
curate model so long as the quality factor of each oscil- 
lator is much greater than unity. For a single oscillator 
with annihilation operator c, this master equation is |24j 

p=-(r/2)(n T + l)({ C t c ,p} + -2 C p C t) 

-(r/2)n T ({cc\p} + ~2Jpc), (4) 

where p is the density matrix, T is the damping rate, 
and c is the annihilation operator for the oscillator. The 
constant is the average occupation number of the os- 
cillator at the ambient temperature, T, and is given by 
riT = \/[exp{—hv/kT) — 1], where k is Boltzmann's con- 
stant and v is the angular frequency of the oscillator. 
Both oscillators are subject to this master equation, with 
v and r replaced with the appropriate quantities. 
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FIG. 1. (a) A simulation of cooling a (target) nano-resonator 
with frequency u), using a frequency-conversion (modulated) 
coupling to an (auxiliary) L-C oscillator. The initial thermal 
state of the target has occupation (n) = 3.68. The solid lines 
have coupling g = uj/(20n) = 2k, where k is the damping 
rate of the auxiliary. The Q values are 10 4 (light grey), 10 5 
(medium grey), and 5 x 10 5 (black). The dashed lines have 
the same damping rate for the target as the black solid line, 
but now with oj/g = 407r (grey), and uj/g — > oo (black). The 
straight dashed lines are the simple steady-state estimates 
provided by Eq.|l]) for the three different damping rates, 
(b) A plot showing the cooling achieved by a single state- 
swap. The parameters for the system are those used for the 
medium-grey solid line in Fig. [T^,. The dashed-line, for which 
the initial occupation number is (n) = 3.68, is merely a blow- 
up of the medium-grey solid line in (b). The solid line has 
(n) = 20. 

We now simulate the evolution of the target and auxil- 
iary, under the interaction given by Eq.([3]), and the ther- 
mal damping of both oscillators. To perform this simula- 
tion we use the Monte Carlo wave-function method [21] , 
combined with a Milstein stochastic integrator We 
find that 4096 trajectories are sufficient to obtain accu- 
rate results. We are somewhat restricted in the size of 
the Fock space we can use for the simulation, but this is 
not especially important since it is the cooling factor that 
we wish to confirm. Previous results on side-band cool- 
ing, consistent with the simple expression derived above, 
confirm that the cooling factor is largely independent of 
the ambient temperature. 

For our first set of simulations we choose an initial 
thermal state for the target with occupation number 
(n) = 3.68. This corresponds, for example, to a 100 
MHz resonator at 20 mK, and can be simulated with 32 
Fock-states for each resonator. In Fig. [T^, (solid lines) 
we plot the evolution of (n) for the target for a 20 MHz 
resonator with g = 2k = 2 MHz, and three values of 
the target damping rate: r y/2ir = 2 KHz, 200 Hz, and 
40 Hz (corresponding to quality factors Q = 10 4 , 10 5 , 
and 5 x 10 5 ). We also plot the steady-state value for (n) c 



as predicted by the simple formula in Eq.([T]) (straight 
dashed lines) . Here we see that the frequency-conversion 
interaction swaps the phonons between the nanomechan- 
ical resonator (the target) and the auxiliary oscillator. 
The auxiliary's large damping rate cools the system un- 
til the swapping ceases. We see that for Q = 10 4 the 
steady-state value is within 30% of the simple estimate. 
The efficacy of the estimate appears to degrade as Q is 
reduced, but this is not the case — as we increase Q, and 
thus the cooling factor, we are instead seeing the limit- 
ing effect of the size of g/ui (the failure of the rotating 
wave-approximation). To show this, we also plot (n(t)} 
with the time-dependent part of the interaction dropped 
(the black-dashed curve), and we see that the resulting 
steady-state is once again within 30% of that predicted 
by Eq. ([!]) . This shows us that our simple formula consis- 
tently overestimates the cooling by a factor of « 1.3, so 
long as g/uj is sufficiently small. From these plots we also 
see that the finite value of g/u does not impose a fixed 
limit on the achievable temperature, but as has already 
been noted in previous analyses, acts like an additional 
heating rate JTJ [T2] . 

We now examine the state-swapping induced by fre- 
quency conversion, which is already evident in Fig. [l^i. 
Our numerical simulations show that the damping of the 
auxiliary, k, has little effect on the fidelity of the swap, 
even when it is the same order-of-magnitude as the inter- 
action rate. It does, however, shift the time at which the 
swap occurs. For example, for k = g/2 the swap time is 
1.1777r/(2<7). Note that while cooling in the steady state 
requires k 7, ground-state preparation using a state- 
swap does not. Thus if the ratio ft/7 is low enough, the 
state-swap will achieve a lower temperature than side- 
band cooling. The degree to which the state-swapping 
can prepare the resonator in the ground state (that is, 
how cold it can cool) in this case depends upon the effi- 
cacy of the rotating wave approximation (g <C w) and the 
heating of the target that takes place during the swap. 
The latter is set by the ratio of g to the heating rate 771^- 

In Fig. [TJd we show the cooling achieved by a state- 
swap for g — lu/(20it) = 2k, Q — 10 5 , and two values 
for the temperature of the target resonator, nx = 3.68 
and 20. For these values of the temperature we obtain 
cooling factors of 725 and 758, respectively, showing that 
the cooling factor varies only slowly with the ambient 
temperature. By increasing the ratios w/g and 5/7 we 
can obtain higher cooling factors. An example of this is 
shown by the dashed-lines in Fig. [l^, in which the first 
minimum of (n) is less than 10^ 3 . 

Coherent feedback control: To control the coherent 
state of the target — that is, control its phase and ampli- 
tude — first consider when the target is at zero temper- 
ature. In this case, driving the auxiliary oscillator with a 
classical input signal will in turn provide a coherent drive 
for the target via the coupling with the auxiliary. The 
resulting steady-state phase and amplitude of the target 
is thus determined by the drive on the auxiliary. Since 
the thermal noise affecting the target is uncorrelated with 
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the drive, the resulting steady-state of the target is sim- 
ply the cooled thermal state, displaced in phase space to 
the amplitude and phase dictated by the drive. The aux- 
iliary thus simultaneously cools and controls the target. 

The coherent control is easily analyzed by using the 
Heisenberg equations of motion for the coupled oscilla- 
tors. To do this we make the RWA, because while the mo- 
tion can be exactly solved if we keep the time-dependent 
terms in the Hamiltonian (Eq.([3|), the resulting expres- 
sions are very complex. The Heisenberg equations for the 
coupled oscillators, in the RWA are 

±( a \ _ _( i ig \ ( a\ ( Train \ (e s 
dt\b)- \ig %)\b) + \ s/K[b in + P] ) (5J 

Here /? is complex, and gives the phase and amplitude 
of the signal driving the auxiliary resonator. The mag- 
nitude of j3 is related to the mi crowave power driving 
the auxiliary by |/?| = y / P/(hQ) [26] ■ The symbols Oi n 
and &i n are mutually independent quantum white-noise 
sources, characterized by the correlation functions |26j 

(a in (t)al(t + T)) = (l + n T )S(r) (6) 
(al(t)a in (t + T))=n T 5(r), (7) 

and (bin{t)bl n (t + r)) = <5(r). Solving these equation for 
a(t) and b(t), and setting 7 <C k and g < k/4, the re- 
sulting steady-state coherent amplitude of the target res- 
onator is 




Thus the phase and amplitude of the target is controlled 
by the auxiliary. The rate at which the target responds to 



changes in the driving is gi ven by the two decay constants 
A± « k/2 ± ^k 2 /4-V ~ k. 

The steady-state phonon number is |(a) ss | 2 , plus the 
contribution (n) c due to the temperature set by the cool- 
ing. While (n) c is approximately given by our simple ex- 
pression in Eq.([T]), the more exact expression, even within 
the RWA is rather complex: 

where N± = (A| + g 2 ) 1/2 . 

To summarize, we have shown that a frequency- 
converting, linear coupling between a high and low- 
frequency resonator can be used to cool the latter to its 
ground state using a state-swap. Further, in the steady- 
state, this process of state-swapping becomes a coherent 
feedback control loop, in which the mechanism of noise 
reduction corresponds to that of sideband cooling. This 
connection between state-swapping and coherent feed- 
back control may provide insights into coherent control 
of more complex systems. 
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